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ABSTRACT 

We develop the pseudo-C^ method for reconstructing the Cosmic Microwave Back- 
ground (CMB) temperature and polarization auto- and cross-power spectra, and estimate 
the pseudo-CV covariance matrix for a realistic experiment on the cut sky. We calculate 
the full coupling equations for all six possible CMB power spectra, relating the observed 
pseudo-CVs to the underlying all-sky C/s, and test the reconstruction on both full-sky 
and cut-sky simulated CMB data sets. In particular we consider the reconstruction of 
the Ci from upcoming ground-based polarization experiments covering areas of a few 
hundred deg 2 and find that the method is fast, unbiased and performs well over a wide 
range of multipoles from I — 2 to £ — 2500. We then calculate the full covariance matrix 
between the modes of the pseudo-temperature and polarization power spectra, assuming 
that the underlying CMB fields are Gaussian randomly distributed. The complexity of 
the covariance matrix prohibits its rapid calculation, required for parameter estimation. 
Hence we present an approximation for the covariance matrix in terms of convolutions 
of the underlying power spectra. The coupling matrices in these expressions can be es- 
timated by fitting to numerical simulations, circumventing direct and slow calculation, 
and further, inaccurate analytic approximations. We show that these coupling matrices 
are mostly independent of cosmology, and that the full covariance matrix for all six 
pseudo-C^ power spectra can be quickly and accurately calculated for any given cosmo- 
logical model using this method. We compare these semi-analytic covariance matrices 
against simulations and find good agreement, the accuracy of which depends mainly on 
survey area and the range of cosmological parameters considered. 

Key words: methods: data analysis - methods: statistical - cosmic microwave back- 
ground - large-scale structure of Universe 



1 INTRODUCTION 

The study of the CMB has undergone a revolution in the past 
few years, culminating in the first-year Wilkinson Microwave 
Anisotropy Probe (WMAP) observations (Bennet et al. 2003). 
The WMAP results along with those from large scale galaxy 
surveys such as the 2-degree Field Galaxy Redshift Survey (2dF- 
GRS; Percival et al. 2001) and the Sloan Digital Sky Survey 
(SDSS; Tegmark et al. 2004) have revolutionised our under- 
standing of the Universe, yielding measurements of cosmologi- 
cal parameters to unprecedented accuracy and providing robust 
tests of the standard cosmological model. Furthermore, the pace 
of progress in observations of the CMB continues unabated. 
Following on from the exquisite WMAP observations of the 
CMB temperature field, there is a number of ground-based and 
balloon-borne experiments designed to measure the polarization 
of the CMB radiation. Indeed the even-parity E-mode polariza- 
tion signal has already been detected, initially by the Degree 
Angular Scale Interferometer (DASI) experiment (Kovac et al. 
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2002). The first-year WMAP observations then yielded a detec- 
tion of the cross correlation between the temperature and the 
E-mode polarization (Kogut et al. 2003) and very recently, the 
Cosmic Background Imager (CBI) has detected both the E-mode 
signal and its cross correlation with the temperature field (Read- 
head et al. 2004). 

Over the next few years, one of the central goals of CMB 
research will be to tie down the polarization power spectra and 
take advantage of their different cosmological dependencies, 
particularly their use as a direct probe of inflation (e.g. Dodel- 
son, Kinney & Kolb, 1997; Kamionkowski & Kosowsky, 1998). 
With this in mind, a precise measurement of the CMB polar- 
ization field would be extremely useful and there are a number 
of upcoming experiments designed to make such measurements 
of both the E-mode and, perhaps more importantly, the odd- 
parity B-mode polarization signal. In addition to planned satel- 
lite missions with the capability of large-scale, high-precision 
polarization measurements (e.g. the Planck, Tauber 2004, and 
proposed Inflation Probe 1 experiments), there are a number of 

1 http://universe.gsfc.nasa.gov/program/inflation.html 
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funded ground-based experiments designed to measure the po- 
larization field with very high signal-to-noise over a limited re- 
gion of sky. In particular, the QUaD 2 (QUEST 3 and DASI), and 
BICEP 4 (Background Imaging of Cosmic Extragalactic Polar- 
ization) experiments should begin observations in 2005, while 
the CLOVER (CMB Polarization Observer; Taylor et al. 2004) 
experiment is planned for 2008. A detection and measurement 
of B-mode polarization in the CMB would provide a major step 
forward in the field and, if a B-mode signal were found to be 
present on large enough angular scales, this would constitute 
a direct detection of primordial gravitational waves (Seljak & 
Zaldarriaga 1997; Kamionkowski, Kosowsky & Stebbins 1997). 
Such a measurement would have immense consequences for the- 
ories of the early Universe, possibly providing a handle on the 
energy scale of inflation (e.g. Kinney 1998). 

These ground-based experiments will measure the CMB 
signal over finite areas of sky, typically a few hundred square 
degrees. There already exist well-established methods for re- 
covering the true underlying CMB temperature power spectrum 
from observations of such finite areas of sky and these are well 
tested (Gorski 1994; Tegmark 1997; Bond, Jaffe & Knox 1998; 
Dore, Knox & Peel 2001; Szapudi, Prunet & Colombi 2001a; 
Szapudi et al. 2001b; Wandelt, Hivon & Gorksi 2001; Hansen, 
Gorski & Hivon 2002; Hivon et al. 2002). Although these meth- 
ods, in principle, should be equally applicable to the polarization 
power spectra, many have not been tested with detailed sim- 
ulations. Hansen & Gorski (2003) investigate temperature and 
polarization power spectra estimation using Gabor transforms 
while Chon et al. (2004) have extended the analysis of Szapudi 
et al. (2001a) to construct a polarization power spectrum estima- 
tor based on correlation functions. 

In this paper, we extend, and test with simulations of the 
CMB temperature and polarization fields, the pseudo-C^ method 
for fast estimation of CMB power spectra. Although not opti- 
mal, the pseudo-C^ method has the advantage that it is fast and 
completely applicable to mega-pixel data sets such as those ex- 
pected from upcoming satellite and ground-based experiments. 
We have generalised the method first introduced for the tem- 
perature field, T, by Hivon et al. (2002), and extended to the 
E and B modes of polarization by Kogut et al. (2003), to the 
full set of six CMB power spectra: Cj T , Cf E , Cf B , Cj E ', 
Cj B and Cf B . Note that this method for recovering power 
spectra is completely applicable to other cosmological data sets 
such as large-scale weak lensing surveys, via the transforma- 
tion, {T, Q, U, E, B} — > {/i, 71, 72, k, (3}, where /i is the mag- 
nification, 71 and 72 are the two components of the observed 
shear, k is the even-parity lensing convergence and j3 is a non- 
gravitational odd-parity distortion. 

We also present, and test with simulations, a new method 
for fast and accurate calculation, for any given cosmological 
model, of the full covariance matrix of all six pseudo-C^ power 
spectra. This method will be extremely useful in estimating 
cosmological parameters directly from the observed pseudo-CV 
power spectra. We also present the analytic expressions for the 
covariances (see also the calculations of Hansen & Gorski 2003 
and Challinor & Chon 2004). These expressions are not easily 
simplified and evaluated so we propose fitting the mixing matri- 
ces which appear in these expressions from Monte-Carlo (MC) 
simulations. Provided that these mixing matrices are mostly in- 
dependent of cosmology, this method allows one to calculate the 

2 http://www.astro.cf.ac.uk/groups/instrumentation/projects/ 
1 Q and U Extragalactic Submillimetre Telescope 
4 http://bicep.caltech.edu 



covariance matrix with relative ease for any given cosmological 
model. 

The paper is organised as follows. In Section 2, we present 
the calculation of the pseudo-CV method relating the observed to 
full-sky power and cross-power spectra for the temperature (T) 
and polarization (E, B) fields of the CMB. In Section 3, we cal- 
culate the full covariance matrix of the pseudo-CV estimator and 
we outline our method for calculating it for a given cosmolog- 
ical model. In Section 4, we test the pseudo-C^ reconstructions 
and the covariance matrix fitting method on realistic simulations 
of CMB temperature and polarization experiments for both full- 
sky and finite-area data sets, and we also discuss the results and 
implications of these simulation tests. Our conclusions are pre- 
sented in Section 5. 

2 PSEUDO C e METHOD FOR CMB POWER SPECTRA 

As is well known, in the currently favoured inflationary based 
cosmological model, the density perturbations are expected to 
be isotropic and nearly Gaussian distributed (e.g. Liddle & Lyth 
2000). In this case, all the information present in the CMB tem- 
perature and polarization fields can be compressed without loss 
into the power spectra of the fields. The extraction of these 
power spectra from pixelized and cut-sky maps has become a 
real challenge when analysing large CMB data sets. Indeed, the 
standard approach, based on maximum likelihood methods, ei- 
ther in Fourier sapce (Gorski 1994; Gorski et al. 1994) or real 
space (Bunn 1995; Bond et al. 1998), is extremely CPU inten- 
sive for a large number of pixels. The need for faster methods 
has led to the development of a number of alternative techniques 
for power spectrum estimation. Szapudi et al. (2001) and Chon 
et al. (2004) have developed power spectrum estimators based 
on correlation functions while Dore et al. (2001) have proposed 
a maximum likelihood estimator based on hierarchical decom- 
position of a CMB map into sub-maps with varying degrees of 
resolution. Hansen et al. (2002) and Hansen & Gorski (2003) 
use Gabor transforms to recover power spectra from apodized 
regions of the sky, while recently, Wandelt, Larson & Laksh- 
minarayanan (2004) have proposed an exact method for recov- 
ery of power spectra from CMB data using Gibbs sampling. An 
overview of power spectrum estimators for the CMB tempera- 
ture field can be found in Efstathiou (2004). 

The application of the pseudo-C^ method to CMB observa- 
tions was first proposed by Wandelt et al. (2001) and Hivon et al. 
(2002). This technique, which is based on the original method 
of Peebles (1973), relates the observed power spectra, Ce, to 
the true underlying spectra, Ce, via a coupling matrix which de- 
scribes the effect of the sky-cut applied to the data. In the follow- 
ing section we begin by deriving analytical expressions for the 
full set of two-point self-correlations and cross-correlations be- 
tween the temperature, T, and polarization, E and B, on a finite 
area of the sky. As we expect theoretically the cross-correlations 
between different parity fields to be zero (B has the opposite 
parity of T and E), any non zero measurement of Cj B or Cf B 
would be indicative of residual noise left in the data, providing a 
powerful test of systematic effects in the experimental pipeline. 
Alternatively, if one is confident that systematic effects are small 
with respect to any measurement of Cj B or Cf B , then these 
spectra can be used to test for parity violations (e.g. Lepora 
1998; Lue, Wang & Kamionkowski, 1999). 

2.1 The CMB temperature and polarization fields 

The CMB radiation is completely characterised by its tempera- 
ture anisotropy, T, and polarization, P, in each direction in the 
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sky. The temperature may be conveniently expanded in spherical 
harmonics, 



T(fi) = ^T em Y lm (Sl) 



lm 



with the inverse transform, 



T lm = dnT(n)Y; m (n), 



(i) 



(2) 



where the integration is taken over the whole sky. As T is a real 
field, the harmonic modes obey the hermicity relation, T/ m = 

(-l) m T,,_ m . 

The polarization field can be described by the Stokes pa- 
rameters, Q and U, with respect to a particular choice of a co- 
ordinate system on the sky (in relation to which the linear po- 
larization is defined). One can conveniently combine the Stokes 
parameters into a single complex quantity representing the polar- 
ization, P = Q+iU. Due to its rotation properties, P — > Pe l2< ^, 
the polarization is a spin-2 field (see Zaldarriaga & Seljak 1997). 
One may thus expand P and its complex conjugate in terms of 
the spin- weighted spherical harmonics, s Yg, m (where s denotes 
the spin), as 

P(Q) = Q(Q) + iU(Q) 

lm 



P*(Q) = Q(Sl)-iU(Sl) 



-2Ye m (Jl), 



(3) 



where the summation in m is over —I < m < I. The E lm 
and Be m are the spin-2 harmonic modes of the so-called electric 
(i.e. 'gradient') and magnetic (i.e. 'curl') components of the po- 
larization field respectively. Using the orthogonality of the spin-s 
spherical harmonics over the whole sphere, 



/ 



dQ s Yt m (Q) 3 Yg* m /(f2) = <5«/<5„ 



(4) 



where the spin states must be equal, the harmonic modes of the 
E and B fields can be found directly from the polarisation field, 
P\ 



E lm = 1 JdQ [p(si) 2 Y; m {o) + p*(rt) 



(n)] 



Bt m = - 



dn [p(si) 2 Y e * m (n) - p*(n) -jr; m (si)]. (5) 



2.2 Finite Area Fields 

In this paper, we are interested in the effects of limited sky cov- 
erage on the measurement of the auto- and cross-power spectra 
of the CMB temperature and polarisation fields. In the case of 
the temperature field, the effect of a finite window function is 
multiplicative, 

f (n) = w T (n)T(n), (6) 

where Wt(0,) defines a weighting function over the whole sky, 
between and 1, and T is the temperature field on the cut sky. 
Outside of the survey, Wt = 0. The actual value of the window 
will depend on the survey scan strategy, foregrounds, etc. The 
effect of the survey window function on the harmonic modes is 
a convolution, 



(7) 



where we have defined the harmonic space window function, 

W%r' = JdSloY tlml (fl)W T (fl) Y; m (fl). (8) 

In the case of the polarization field, the effect of the window 
function is 



p(n) = w P (n)P(n), 



(9) 



where P is the cut-sky field. Note that the temperature and 
polarization window functions will not, in general, be equal; 
Wr(fl) ^ Wp(ft). Using equations (5), we find the cut-sky 
E and B fields are given by (Lewis, Challinor & Turok 2002) 

E lm = y^(-El'm'W«/ mm / + Bl> m 'W u , mm ,) 
I'm' 

Bim = ^(.ZWW+ mm , - E e , m ,W e - e , mm ,), (10) 

I'm' 

where we have defined 

1 . 



Wu'mm' - - 2 ^r —2W%r), (11) 

with the spin-weighted harmonic space window functions (for 
s = ±2) given by 



WP t r = dn 3 Y e/m ,(Q)Wp(Q) s Y; m (Q). 



(12) 



Thus, the effect of observing the CMB radiation over a limited 
area of the sky induces a contamination of the electric signal by 
the magnetic component and vice- versa. This was to be expected 
as the electric/magnetic (or 'gradient/curl') decomposition of the 
polarization vector on a finite patch is non-unique (see Lewis et 
al. 2002). 

We can simplify these expressions for the temperature and 
polarization fields by defining the vector field 



8<m — (Q?m: ©lm) ©lm) — (Te m , E ( m , B £ m ) . 



(13) 



The cut-sky vector, ®e m , is then related to the underlying 
whole-sky field, ®t m , by 



where the window matrix is given by 

/ W%, mml 



tl' mm' 






~^ee'mm' ^U'mm' 



(14) 



(15) 



with W£ w =oW7? m ' given by equation (8). 



2.3 Statistical properties of the CMB on the cut sky 

One can use equations (14) and (15) for the pseudo-temperature 
and polarization fields to relate the values of the six pseudo- 
power spectra of the cut-sky fields to their full-sky counterparts. 
Taking the expectation values over the sphere of the harmonic 
coefficients of the full-sky vector field, ®tm, we obtain the full- 
sky power spectra which we can write as a matrix relation, 



2£ + 1 ^ 



<<=w©L), 



(16) 



where 
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({Ti m Tf m ) (Tt m E* tm ) {T lm Bl m )\ 
(E em Tt m ) {E tm E% m ) (E lm B* lm ) J . (17) 
{Bt m Tt m } (Bt m E* tm ) (B im B* lm )J 

For fields defined only over finite areas of the sky, we can use 
equation (14) to obtain an expression coupling the cut-sky power 
spectra to the full-sky spectra, 



21 



m 

TT £ E w^'c^w^Y, (18) 



21 

V mm' 

which, for clarity, can be written as 

C e — M«/C^/. 



(19) 



In Appendix A, we outline the derivation of the full coupling 
matrix, M. Armed with these relations, the true auto- and cross- 
power spectra can be estimated from the pseudo-power spectra 
by inverting equation (19), 



(20) 



In Section 4, we shall apply this relation to a set of numerical 
simulations of the CMB. But before we proceed, we now turn 
our attention to calculating the covariances between the pseudo- 
C't spectra. 

3 COVARIANCE MATRIX OF THE PSEUDO-C £ 
ESTIMATORS 

In addition to estimating the full-sky temperature and polariza- 
tion power spectra from the cut-sky spectra, we would also like 
to understand the covariances between the different pseudo-CV 's 
for model fitting or parameter estimation. In the case of maxi- 
mum likelihood parameter estimation, the covariance matrix is 
required at each point in parameter space sampled by the param- 
eter estimation procedure. Calculating it from MC simulations 
is then clearly unfeasible if one wishes to investigate a large 
number of parameters. We therefore require an accurate and fast 
method for calculating the covariance matrix of the pseudo-Cf 
measurements. 

The covariance matrix of the various pseudo-C^ spectra can 
be expressed exactly in terms of the components of the window 
matrix, equation (15). Defining the covariance of two measur- 
able quantities, X and Y, as 



(AXAY) = {(X-{X))(Y-{Y))), 



(21) 



where the angled brackets denote an ensemble average, the co- 
variance of the pseudo-Cj? spectra, for Gaussian-distributed tem- 
perature and polarization fields, is given by (c.f. Hansen & 
Gorski 2003; Challinor & Chon 2004) 

(AC? Y AC^ N ) = 



E( r ,AD n BC y{XA 



ND,MC,YB} 



n AC n BD V {XA,MC,ND,YB}\ 

where the X«/^ 2 terms are given by 



(22) 



{XA,ND,MC,YB} 



we 1 e 2 



(2£+l)(2£' + l) 

EixrXA ff\rND \*TxrMC l\%r YB \* 

Wu 1 mm 1 (Wt'e 1 m'm 1 ) W V l 2 m' m 2 ( W U 2 mm 2 ) 



(23) 



Here, the final sum is over all m, m! , mi and mi and a sum over 
all possible values of A, B, C and D is implied. The W^?^ m i 
matrices are just the entries of the window matrix, equation (15). 
Explicitly, these are given by (suppressing the subscripts for 
clarity): W TT = W°, W EE = W+, W BB = W+, W EB = 
W~ and W BE = —W~. For all other combinations of X and 
Y, W XY = 0. Similar expressions to equations (22) and (23) 
also arise in the calculation of the covariance matrix of cross- 
power spectra obtained from different channels of the same ex- 
periment for the temperature-only case (Tristram et al. 2005). 
For the temperature-only case, equation (22) reduces to 



(ACf T ACji 1 



TT j^{TT ,TT ,TT ,TT) 



(24) 



In the case of a full-sky data set with a narrow galactic cut ap- 
plied, Efstathiou (2004) has proposed approximating this exact 
expression by replacing Cf^ and Cj 2 with Cj T and Cj, T 
and applying the completeness relation for spherical harmonics 
(Varshalovich et al. 1988). This yields an expression for the co- 
variance matrix of the pseudo-C^ measurements in terms of the 
temperature coupling matrix of equation (A12). Unfortunately 
an equivalent approach for the full set of pseudo-CV estimators 
does not result in a simple expression for the covariance matrix 
since the completeness relation for spin spherical harmonics is 
only valid for spherical harmonics of the same spin. Recently, 
Challinor & Chon (2004) have performed additional approxi- 
mations yielding analytic expressions for a restricted part of the 
pseudo-CV covariance matrix in terms of the Wigner-3j symbols. 

Here we propose a different method for calculating the co- 
variances of equation (22). Generalising the approach of Efs- 
tathiou (2004), we make the symmetrized approximation, 



n AD n BC 



/r'ADr^ADr'BCr'BC 

V £ U f U f ' 



(25) 



(and similarly for the second term in equation 22) and then fit 
the resulting Ce -independent terms to the covariance measured 
from MC simulations. That is, after applying the above approx- 
imation we are left with an equation of the form (where, again, 
summation over values of A, B, C and D is implied), 

{ACf Y ACf N ) « 

/r<ADr'ADr'BCr , BC V {XA,ND,MC,YB} 



Jr^ACrACrBDr'BD ^-{xa,mc,nd,yb} 
where the C^-independent mixing terms are given by 



(26) 



Y {XA,JVD,AfC,yB} _ \ ^ V {XA,ND,MC,YB} 

For the covariances of the temperature pseudo-C^ 's only, we find 

(ACj T ACj, T ) = 2Cj T Cj, T X. { ^, T ' TT ' TT ' TT} , (28) 

and in the case where the weighting function, Wt((2), is a sim- 
ple mask with values of 1 or 0, the X«/ matrix simplifies to 
(Efstathiou 2004): 



X 



{TT,TT,TT,TT} 



1 



21' + 1 



TT,TT 



(29) 



TT TT 

where M u , ' is the temperature coupling matrix of equa- 
tion (A12). Unfortunately, the remaining 20 terms of the X- 
matrix do not simplify so easily. We provide the explicit expres- 
sions for the mixing matrices and the 21 possible covariances in 
Appendix B. Note that the mixing matrices exhibit some sym- 
metries, in paticular, X^, cd = Xf|, 6a . While we could proceed 
to try and find approximations to the X-matrix elements, this 
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would introduce further uncertainties to the estimation of the co- 
variances. Instead, we propose fitting the X-matrix terms to MC 
simulations of the survey under consideration. Our assumption 
is that the mixing terms, X«/, are very nearly independent of 
cosmology and predominantly dependent on the sky-cut applied 
to the data. Note that this same assumption is also required if we 
were to try and find analytic approximations to the X«/ terms. 
We test this assumption and our proposed fitting method with 
simulations in Section 4.2. 



4 TESTS ON SIMULATIONS 

Equations (20) and (A12) - (A 17) provide us with a method for 
recovering the various underlying full-sky CMB power spectra, 
Ci, from a set of pseudo-power spectra, Ci, measured from a 
cut-sky map. In addition, the method outlined in Section 3 pro- 
vides us with an approximation to the pseudo-CV covariances for 
any given cosmological model. In this section, we test these pro- 
cedures on simulated CMB temperature and polarization maps. 
We first turn to the Ct —> Ci reconstruction, after which we 
investigate the covariance matrix approximation. 

4.1 Power spectra reconstructions 

4.1.1 Simulated CMB maps 

We have used the HEALPix 5 package (Gorski, Hivon & Wan- 
delt 1999), synfast to create simulations of the temperature (T) 
and polarization (Q,U) fields of the CMB sky. We have con- 
sidered two experiments: (1) the Planck experiment consisting 
of two sky surveys over 14 months and (2) a combination of 
the QUaD polarization experiment (Church et al. 2003) with 
the 4-year temperature data from the WMAP satellite. The pix- 
elization scheme used in all the simulations is characterised by 
the HEALPix parameter, N^c, where we have used a value of 
A?aide = 2048, giving an average pixel area of 2.95 arcmin 2 . 
Using such a small pixel size presents one with some comput- 
ing speed and data storage issues, especially when it comes to 
running large numbers of MC simulations. However, this value 
of A^ide is necessary to ensure that both the Planck and QUaD 
beams, with full width at half maximums (FWHM) of 5 and 4.2 
arcmin respectively, are adequately sampled. If this was not the 
case, then the resulting temperature and polarization maps would 
not be strictly bandwidth-limited and the ensuing calculation of 
the pseudo-C^ 's from the maps would potentially be biased at 
high I. Uncorrected Gaussian random noise has been added to 
the simulations in order to mimic the instrumental noise effects 
in the detectors. The sensitivity levels used are those quoted on 
the Planck 6 and WMAP 7 websites and, for the QUaD experi- 
ment, in Bowden et al. (2004). In each case, we have combined 
sensitivity levels from the different channels in the three experi- 
ments. The experimental parameters we have used to create the 
simulations are given in Table 1. 

The simulations are random Gaussian realisations of theo- 
retical CMBFAST 8 (Seljak & Zaldarriaga, 1996) temperature and 
polarization power spectra. Note that we have used a slightly 
modified version of CMBFAST in which the amplitude of the ini- 
tial curvature perturbations is characterised by the parameter, A. 

5 available from http://www.eso.org/science/healpix 

6 http://www.rssd.esa.int/index.php ?project=PLANCK 

7 http://map.gsfc.nasa.gov 

8 available from http://www.cmbfast.org 



This is the same parametrization used by the WMAP team (see 
Verde et al. 2003 for details). The cosmological parameters used 
to generate the CMBFAST spectra were the best-fitting parameters 
of the WMAP experiment for their power law ACDM model 
(Spergel et al. 2003): 

{fi m , n A , Ob, h, n s ,T, A} 

= {0.27, 0.73, 0.046, 0.72, 0.99, 0.17, 0.90}. (30) 

Here, Q, m is the total matter density, Qa is the dark energy den- 
sity, n b is the energy density in baryons, h is the Hubble con- 
stant in units of 100 kms _1 Mpc _1 , n s is the spectral index of 
the initial perturbations, r is the optical depth to the last scatter- 
ing surface and A is the parameter used by the WMAP team to 
characterise the amplitude of the initial perturbations. 

In addition, to test the reconstruction of the BB power spec- 
trum, we have included tensor perturbations in our cosmological 
model. We have used a tensor-to-scalar ratio of 

T / 5 = #7FT = a05 < (31) 

where A^(fco) is the power spectrum of the initial curvature 
perturbations and A T (fcn) is the power spectrum of gravitational 
waves from inflation - both evaluated at fco = 0.05 Mpc^ 1 . We 
have also included the effects of gravitational lensing in our cos- 
mological model which increases the amplitude of the B-mode 
polarization signal on small to medium scales. Multipoles from 

1 — to i — 3000 were used in the map-making process. After 
generating the maps, they were smoothed with a Gaussian beam 
with a FWHM of 13 (WMAP), 5 (Planck) and 4.2 (QUaD) ar- 
cmin. These FWHM values correspond to the highest resolution 
channel for each of the three experiments (see Table 1). Exam- 
ple maps of the temperature (Planck) and polarization (QUaD) 
fields are shown in Fig. 1. 

4. 1.2 Reconstructing the power spectra 

To measure the pseudo-CV spectra, we use the HEALPix package, 
anafast to estimate 

6 ' = 2<TiE ito,§ ^ (32) 

m 

where & em = (T<> m , E tm , Bi m ) are the observed spherical har- 
monic coefficients for the temperature and polarization fields. 
These have been calculated from the simulated maps using equa- 
tions (2) and (5) with T(fi) and P(fi) replaced with f(Q) 
and P(fi) respectively. The integrals in these equations have 
been approximated as sums over discrete pixels. For the tests 
we present here, we have used the simplest possible weight- 
ing scheme: Wt = 1 for observed pixels and Wt = oth- 
erwise (and similarly for Wp). Note that the method of Section 

2 also allows for more complicated window functions, e.g. in- 
verse noise weighting for each pixel. The anafast package ex- 
ploits various symmetries of spherical harmonics to reduce the 
computation of the &tm coefficients from a JVpix^max process 
to one which scales as -Wp-^max where 7V p i x is the number of 
pixels in the map and i? max is the maximum multipole to be mea- 
sured. 

When it comes to applying the method of Section 2 to re- 
cover the CMB spectra from the pseudo-C^ 's measured from an 
experiment with a large sky-cut, it turns out that equations (20) 
and (A12) - (A17) cannot be applied directly since the coupling 
matrix, M become singular and un-invertible for large sky-cuts. 
In this case, it is necessary to bin the power spectrum estimates 
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Table 1. Experimental parameters used to create the simulated Planck, WMAP and QUaD datasets. Note that, for the Planck simulations, we have 
considered only the 4 lowest frequency channels of the Planck High Frequency Instrument (HFI) and, for the WMAP simulations we have included the 
Q, V and W bands only. The other bands of these experiments will mainly be used to constrain foregrounds. Noise levels (AT/T) are given separately 
for temperature and polarisation in terms of the sensitivity (lcr) to intensity (Stokes I) and polarised intensity (Stokes Q and U) fluctuations, relative 
to the average temperature of the CMB (2.73 K). These sensitivity levels are given per HEALPix pixel where the pixelisation scheme we have used is 
characterised by the HEALPix parameter, N s ^ c . For Planck, these noise levels are the performance level goals for 2 sky surveys (14 months). The noise 
levels for WMAP are those expected after 4 years of observations and for QUaD, the noise levels are those expected after 2 years observing. Note that 
WMAP will, in fact measure polarisation, and QUaD will measure temperature. However, in this paper, we are interested in the combination of the 
WMAP temperature and QUaD polarisation measurements and so do not quote these sensitivity levels. 







Planck HFI 






WMAP 




QUaD 


Sky fraction (/ sky ) 




0.826 






0.826 




0.0078 


Frequency (GHz) 


100 


143 217 


353 


41 


61 


94 


100 150 


Beam FWHM (arcmin) 


9.5 


7.1 5.0 


5.0 


30 


21 


13 


6.3 4.2 


AT/T X 10 5 (temperature) 


1.07 


0.91 1.40 


4.27 


2.85 


3.47 


3.28 




AT/T X 10 5 (polarisation) 


2.14 


1.73 2.85 


8.67 








0.36 0.32 


HEALPix parameter, N s \^ e 




2048 






512 




2048 


Average pixel area (sq. arcmin) 




2.95 






47.2 




2.95 




Figure 1. Example simulated maps for the Planck (temperature - left) and QUaD (polarization - right) experiments. A ±10° galactic cut has been excised 
from the simulated Planck maps. The temperature field is shown as the colour scale and the polarization field in the QUaD simulation is represented 
by the vector field. The orientation of each vector indicates the direction of polarization while the length indicates the polarization magnitude. The 
simulated QUaD map measures 17.5 X 17.5 square degrees - similar to the proposed QUaD survey area (Bowden et al. 2004). 



and reconstruct the true CMB spectra in terms of bandpow- 
ers. Following closely the temperature analysis of Hivon et al. 
(2002), we calculate full-sky CMB bandpowers (Pt) as 



(33) 



where the subscript b denotes the bands whose boundaries in £- 

space are defined by < ^g h < Mow*'- Pbl is a Mnnm g 
operator in £-space defined by 



2tt ,(M-1) Jb) > 



0. 



if 2 < e (b) < f < f (b+1) 

11 Z Mow =* 1 *- Mow 

otherwise. 



(34) 



In equation (33), we have also introduced the binned coupling 
matrix, M bi / which is constructed from the coupling matrix of 
Section 2 using 



where Qa> is the reciprocal operator of Pbi 



(35) 



e(e+i) > 
0, 



if 2 < ( {b) < t < f (b+1) 

11 Z Mow * ^ Mow 

otherwise. 



(36) 



In this way, we approximate the true CMB spectra as being 
piece-wise constant where the quantity Pj, = £(£ + 1)C^/2tt 
is taken to be constant within each band. We can also include in 
the coupling matrix a correction for the effects of the beam and 
finite pixel size, 



(37) 



where Bi is the combined window function describing the beam 
smoothing and pixelization of the experiment. To calculate the 
Wigner 3j symbols, which are required for the coupling matri- 
ces of Appendix A, we have used the publicly available SLATEC 
fortran subroutine, DRC3JJ.f. This routine uses the algorithm of 
Schulten & Gordon (1975) which makes use of both forward 
and backward recurrence relations to maintain numerical stabil- 
ity. It is both rapid and accurate, even for large multipole values. 
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In order to calculate the errors and covariance of the result- 
ing bandpower measurements, we run a number (~ 100) of MC 
simulations. Note that we find very little difference in the covari- 
ances measured from 50 and 100 MC simulations and therefore 
conclude that our suite of 100 simulations is sufficient to yield a 
good estimate of the covariances. After measuring the P;, band- 
powers from each of the simulations, their covariance is then 
given by the run-to-run scatter among the simulations: 

(AP 6 AP^) = ((P(, - P t )(p; - p' 6 ))mc, (38) 

where P(, denotes the average of each bandpower over all simu- 
lation runs. Finally, the effects of the noise can be removed from 
the final bandpower estimates by replacing equation (33) with 

P 6 = M w Pvi & - (N<)*c), (39) 
v i 

where (N^) mc is the average pseudo-power spectrum measured 
from a new set of 100 signal-free noisy MC simulations. An 
estimate of the full-sky noise power spectrum is then 

N b = M w Pvi Wmc (40) 

v i 



4.1.3 Results 

We have applied the C e — > Ce reconstruction to the simu- 
lated Planck and QUaD + WMAP data sets in terms of the 
bandpower formulation of the previous section. After applying 
a ±10° galactic cut to the simulated Planck and WMAP maps, 
the pseudo-CV spectra were measured up to a maximum mul- 
tipole of £ max = 3000. Equation (39) was then applied to the 
measured Cf's in order to recover the true underlying full-sky 
power spectra. The reconstruction was carried out up to a max- 
imum multipole of £ m ax = 3000. However, for the very high 
i?-modes, there will be significant contributions to the Ce recon- 
struction from Ci with £ > 3000. We therefore only plot the 
recovered Ce up to a maximum of £ = 2500, thus ensuring that 
contributions from Ce with £ > 3000 are negligible. Note that 
the resulting Ce measurements are also corrected for the effects 
of the Gaussian beam and finite pixel size. 

In Fig. 2, we plot the average recovered spectra from our 
suite of Planck-like simulations. For these reconstructions, we 
have used a linear binning scheme with bin sizes of A£ =100 
for the TT, EE, TE and TB spectra and A£ = 300 for the BB 
and EB spectra. We use the coarser binning scheme for the latter 
two spectra due to the lower significance of the signal for these 
spectra. The error bars plotted are the la errors calculated from 
the diagonal elements of the covariance matrix (equation 38). 
Note that we plot the errors appropriate for a measurement from 
a single experiment. The error in the mean recovered value, av- 
eraged over all the simulation runs, is much smaller. We note 
that the recovery of the TT, EE and TE spectra is unbiased over 
the entire range of multipoles plotted (0 ^ £ < 2500) while the 
BB recovery is also unbiased with a significant detection of BB 
power over a restricted range of l?-space. Comparing the recon- 
structed power spectra with the input CMBFAST models, we find 
reduced-x 2 values of 1.86, 0.83, 0.38 and 0.03 for the TT, TE, 
EE and BB spectra respectively. We suspect that the relatively 
poor reduced-x 2 value of 1.86 for the recovery of Cj T from the 
Planck simulations is due to shot noise from the finite number of 
simulations. This shot noise effect is only seen in the recovered 
Planck Cf T because of the very low temperature noise levels 
we have used in creating the Planck simulations. Finally we note 



that the non-cosmological TB and EB spectra are consistent with 
zero over the entire multipole range. 

The recovered bandpowers for the combined QUaD + 
WMAP simulations are shown in Fig. 3. Here, we plot the TT, 
TE and TB spectra up to £ = 1500 only, as beyond this these 
spectra are unconstrained due to the coarser resolution of the 
WMAP simulations (13 arcmin) compared to that of Planck (5 
arcmin). We have used a linear binning of A£ — 75 for the 
TT, TE and TB spectra, A£ = 100 for the EE spectrum and 
A£ = 300 for the BB and EB spectra. Once again, the error bars 
plotted are the la run-to-run errors calculated from the scatter 
among the simulations (equation 38). We find excellent recovery 
of the TT and TE spectra for < £ < 1300. Beyond I = 1300, 
the effects of the 13 arcmin resolution of WMAP are evident in 
these spectra. The recovery of the EE spectrum is excellent over 
the entire range, ^ £ ^ 2500, while for the BB spectrum, the 
recovery is unbiased with a significant detection of the BB power 
for £ < 1200. These levels of detection for the polarization spec- 
tra agree well with the predicted QUaD detection levels from the 
Fisher matrix analysis of Bowden et al. (2004). Again, the non- 
cosmological TB and EB spectra are consistent with zero over 
the entire multipole range. Comparing the recovered bandpow- 
ers with the input spectra, we find reduced-x 2 values of 0.43, 
0.03, 0.62 and 0.03 for the TT, TE, EE and BB spectra respec- 
tively. 

4.2 Testing the covariance matrix approximation 

In Section 3, we have outlined our method for calculating the full 
covariance matrix of the pseudo-CV estimators. Here, we test our 
procedure against the covariance of the pseudo-CVs measured 
from simulated data sets. Our procedure is to fit the X«/ ma- 
trices of equations (B2) - (B22) from the covariance measured 
from a number of MC simulations. Our goal here is to reproduce 
the covariances of the pseudo-CV's measured directly from the 
maps rather than the covariances of the reconstructed bandpow- 
ers of the previous section. Note that, in some cases, simulations 
from at least three different cosmological models are required to 
find unique solutions for all the mixing matrices. 

4. 2. 1 Fitting procedure 

To generate the fits, we run a number (~ 1000) of simulations 
for each of 38 cosmological models whose input cosmological 
parameters are drawn from within the WMAP 3a region of pa- 
rameter space surrounding the best-fit model for the WMAP data 
(Spergel et al. 2003). For each of these cosmological models, we 
then measure the pseudo-CV covariances from the scatter among 
the pseudo-C^ spectra measured from the simulations and then 
use these covariances and equations (B2) - (B22) to fit for the 
various ~X.ee> mixing matrices, knowing the input model spec- 
tra used. We average these mixing matrices over the 38 cos- 
mological models to produce the final mixing matrices, which 
we then use to calculate our semi-analytic covariances. In order 
to test the performance of the semi-analytic fitting method, we 
conduct new sets of (~ 1000) simulations using different input 
cosmological parameters from those that were used to fit for the 
~X.ee' matrices. We can then compare the covariance estimated 
from the scatter of the pseudo-CV's measured from the simu- 
lations and the covariance calculated using our fitting method. 
We perform these tests for a number of cosmological models 
spanning the WMAP 1, 2 and 3a confidence regions of parame- 
ter space. We do this for both a Planck-type experiment (with a 
±10° galactic cut excised from the simulated maps) and for the 
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Figure 2. The CMB power spectra and cross power-spectra recovered in terms of bandpowers from the simulated Planck data sets. The input CMBFAST 
spectra are shown as the black curves. The light grey curves are the pseudo-CV spectra as measured by the HEALpix package, anafast and the points 
with error bars are the recovered CMB bandpowers. The spectra are, clockwise from top-left: TT, EE, TE, EB, TB and BB. 



proposed survey design for QUaD with a survey area of ~ 300 
deg 2 . Note that due to the large number of simulations required, 
we have considered only multipoles between ^ I ^ 500 as 
estimating the pseudo-CVs up to I — 2500 from such a large 
number of maps is computationally unfeasible at present. 

The CVrange corresponding to the ensemble of cosmolog- 
ical models we have used to perform the fits are shown as the 
shaded regions in Fig. 4. Also indicated on these plots are the 
WMAP default cosmological model and two comparison mod- 
els drawn from the 1 and 2a WMAP parameter space regions 
which we use to assess the performance of our method (see Sec- 
tion 4.2.3 below). 

A potential issue arises in the fitting procedure due to the 
fact that the Cj E power spectrum crosses zero and is negative 
in some regions of ^-space. A consequence of the symmetrized 
approximation, equation (25), is the appearance of pre-factors 
containing the quantity */ Cf E Cj, E in some of the covariance 
expressions (equations B2 - B22). For some combinations of 
£ and £' then, we are left with imaginary quantities for these 
pre-factors. To avoid this problem, we have simply replaced 
^Cj E Cj E by v / \Cf E Cf E \ in equations (B2) - (B22). In 
addition, for the covariances which depend on a single term con- 
taining Cj E ', some of the mixing matrices have discontinuities 
where Cj E is very close to zero, with resulting discontinuities 
in the fitted covariances. For the small fraction of (£,£')-space 
where this occurs, we have calculated the covariance by linearly 



interpolating from neighbouring parts of the fitted covariance 
matrix. 

4.2.2 Including noise 

For the tests we have conducted in this section, we have included 
Gaussian random noise in the simulations with sensitivity levels 
corresponding to those appropriate for the Planck and QUaD ex- 
periments. To a good approximation, the effects of noise can be 
accounted for in our fitting procedure by replacing the model 
power spectra, C* Y of equations (B2) - (B22) with the sum 
of the cosmological plus noise power spectra, Cf" + N^ Y . 
Here, NjF Y are the full-sky power spectra of the noise which 
we have evaluated from noise-only full-sky simulations. Alter- 
natively one could use equation (40) to reconstruct the full-sky 
noise power spectra from an experiment with limited sky cover- 
age. Note that for the uncorrelated noise we have used in these 
tests, Nj E = Nj B = N EB = 0. The power spectra of the 
noise levels simulated are shown as the heavy dashed lines in 
Fig. 4. Note that for both a Planck- and QUaD-like experiment, 
Nj T is negligible for I < 500. 

4.2.3 Results 

To display the results of the covariance comparisons, we use the 
correlation matrix of the pseudo-CV measurements. Inserting the 
pseudo-Cf spectra into the 6 ^ max -element vector, , which in- 
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Figure 3. The CMB power spectra and cross power-spectra recovered in terms of bandpowers from the simulated WMAP + QUaD data sets. The input 
CMBFAST spectra are shown as the black curves. The light grey curves are the pseudo-Cf spectra as measured by the HEALpix package, anafast and 
the points with error bars are the recovered CMB bandpowers. The spectra are, clockwise from top-left: TT, EE, TE, EB, TB and BB. 



eludes the six pseudo-spectra, in order, TT, EE, BB, TE, TB and 
EB, then the 6 £ max X 6 ^ max correlation matrix is defined by, 

corr(C,,Q,) = = , (41) 

v /(AC f AC,)(AC r AC r ) 

which can take values between —1 (fully anti-correlated) and +1 
(fully correlated) where the diagonal elements, corr(Cf, Ct) = 
1 by design. Note that, for clarity, we plot the correlation matrix 
for multipoles in the range 1 $C i ^ 100 only. However, the 
levels of agreement in this range are typical of the agreements 
found up to <max = 500. To investigate further the agreement 
between the true and semi-analytic covariances, we also plot the 
square root of the diagonal elements of the covariance matrix 
(i.e. the lcr uncertainties) for each of the six pseudo-C^ spec- 
tra. We plot these uncertainties up to the maximum multipole, 
^max = 500 in each case. 

Fig. 5 compares the correlation matrix from the Planck-like 
simulations and from the semi-analytic approximation for a sam- 
ple cosmological model from within the WMAP la parameter 
space region. The agreement between the correlation matrices is 
excellent with the semi-analytic approximation reproducing all 
the significant off-diagonal features seen in the covariance mea- 
sured from the simulations. Note that the covariance from the 
simulations also contains shot-noise due to the finite number of 
simulations used which is partly why the semi-analytic approxi- 
mation appears smoother in these plots. This effect, of course, is 



not important and would disappear with an increase in the num- 
ber of MC simulations used. The lcr errors on the pseudo-Ci 
estimates are plotted in Fig. 6. Again, we find excellent agree- 
ment between the simulations and the semi-analytic fits over the 
entire multipole range. Once again, there are shot-noise effects 
in the simulation curves. We have also tested the semi-analytic 
approximation with simulations for cosmological models on the 
edge of the WMAP 2 and 3a parameter space regions. For each 
of these regions, we have compared with 12 different cosmo- 
logical models. As expected, the further in parameter space one 
strays from the models used to generate the fits, the worse the 
semi-analytic approximation becomes. Even so, we found that 
this is quite a gradual effect and the approximation is still per- 
forming well in the 3a region. This is encouraging as it implies 
that the mixing matrices of equations (B2) - (B22) are fairly in- 
sensitive to changes in the cosmological model, at least for a 
modest sky-cut such as the ±10° galactic cut we have used in 
these tests. 

The effect of a more drastic sky-cut can be seen in Figs. 7 
and 8 which show the comparison between simulations and the 
semi-analytic approximation (again for a sample model from 
the WMAP la parameter region) for the QUaD-like simula- 
tions with a survey area of ~ 300 deg 2 . As expected, this sky- 
cut induces much stronger off-diagonal structure in the correla- 
tion matrices. Most of this structure is reproduced in the semi- 
analytic approximation, albeit slightly less accurately than for 
the case of the near full-sky simulations. As in the case of the 
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Figure 4. The range of CMB spectra used to perform the semi-analytic fitting of the covariance matrix. The shaded regions show the Cf range 
corresponding to the WMAP 3cr parameter space region used for the fits. The best-fit WMAP model is shown as the solid curve. Also shown are two 
example comparison models taken from the WMAP 1 (dashed) and 2cr (dotted) regions of parameter space which we use to demonstrate the performance 
of the fitting method (see Section 4.2.3). The heavy dashed lines show the power spectrum of the noise used in the simulations. Note that for the Planck- 
and QUaD-like experiments which we have considered, the TT noise power spectrum is negligible for £ < 500. 



Planck-like simulations, we have also compared with simula- 
tions from the WMAP 2 and 3a regions. Once again, we found 
the agreement gradually worsening as one moved away from the 
WMAP la region. An example of the performance of the ap- 
proximation outside the WMAP la region is given in Figs. 9 and 
10 which show the correlation matrices and la errors respec- 
tively, for one of the WMAP 2a region models, for the QUaD- 
like simulations. The levels of agreement seen in Figs. 9 and 10 
are typical of the performance of the approximation outside of 
the WMAP la region for a QUaD-like survey area, while for the 
Planck-like simulations, the agreement is even better. 

Figs. 11 and 12 demonstrate the performance of the semi- 
analytic approximation for the off-diagonal elements. These 
plots, which are for the WMAP la region model examples 
shown in Figs. 5 to 8 are for the Planck-like (Fig. 1 1) and QUaD- 
like (Fig. 12) simulations respectively. In each case, we plot the 
correlation matrix as estimated from the simulations and from 
the semi-analytic approximation, and we also plot the residu- 
als, corr(MC) — corr(SA) between the two. The root-mean- 
squared (raw) residuals between the simulated and fitted covari- 
ances are 0.033 and 0.038 for the WMAP la region Planck-like 
and QUaD-like simulations respectively. It is clear from these 
plots that the approximation is reproducing all the significant 
features of the pseudo-C"^ covariances. The remaining correla- 
tions apparent in the residual plots are mostly due to the finite 



number of simulations used and can be reduced to below any 
desired level of significance simply by increasing the number of 
MC simulations. 

Approximate formulae, such as the "Knox" formulae 
(Knox 1995; Jungman et al. 1996) have long been used to 
estimate the power spectrum error bars. In such formulae, 
the fractional error on the power spectrum typically scales as 
y/ 2/(21. + l)/sk y , where / sky is the fraction of sky observed. 
These formulae can also account for the effects of instrument 
noise. They do not, however, take account of mode-mode cou- 
pling or mixing between E and B modes due to finite areas of 
sky. It is interesting therefore to compare the errors predicted by 
such formulae to the errors we have measured from the MC sim- 
ulations and the semi-analytic fits. We have used the expressions 
given in equations (4-11) of Eisenstein, Hu & Tegmark (1999) 
to calculate Knox errors for the six power spectra. These ex- 
pressions predict the errors on the true spectra and so, we have 
converted them to errors on the pseudo-C* spectra using 

(AC e f = J2 M u'(^C e ,f, (42) 

v 

where M«/ are the coupling matrices of Section 2. The errors 
predicted using these formulae are compared to the simulated 
and semi-analytic errors in Figs. 6, 8 and 10 for the Planck- 
and QUaD-like simulations. Although, the approximate formu- 
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lae systematically underestimate the errors over most of the £- 
range, for multipoles above £ ~ 100, the discrepancy is more 
or less constant for a given experiment. In the case of Planck, 
the Knox errors under-estimate the true errors by about 10 per 
cent whereas for the QUaD simulations, the errors are under- 
estimated by ~ 25 per cent. 

4.2.4 Implications for parameter estimation 

We have tested our approximation for the covariance matrix over 
a restricted range of cosmological parameter space and have 
found it reproduces all the significant structure of the covariance 
measured from simulations. 

Ultimately, one wishes to use the semi-analytic approxima- 
tion for parameter estimation and there are a number of options 
available. Firstly, parameter estimation could be performed di- 
rectly on the pseudo-C, measurements. In this case the semi- 
analytic covariance matrix could be used as it is and the model 
spectra would need to be convolved with the coupling matrices 
of Section 2. Alternatively, one could perform the estimation us- 
ing the true model spectra in which case, the pseudo-CV covari- 
ance would need to be converted to full-sky Ct covariances. For 
a small-scale experiment like QUaD, one needs to reconstruct 
the full-sky spectra in terms of bandpowers and so the covari- 
ance of the bandpowers would be needed. This is easily achieved 
using 

(AP b AP ft ,) = 

P bl e (AC, AC,,} (P b2l ,) T (M-,l 2 ) T , (43) 

where M. bb i and P w are the binned coupling matrix and the bin- 
ning matrix respectively, defined in section 4.1.2. 

An estimate of the uncertainties in parameter estimates re- 
sulting from our approximation is clearly also required. We are 
currently investigating with simulations what these uncertainties 
are and our conclusions regarding the effect of the covariance 
approximation on parameter constraints will be presented in a 
future paper. However, it is clear already that the accuracy of 
our approximation is a function of a number of factors. Firstly, 
as indicated by the difference in the two sets of simulations 
we have conducted, the performance of the approximation im- 
proves as the sky coverage of the data increases. This is to be 
expected, of course, since as we move closer to full-sky cover- 
age, the approximation converges on the theoretical covariance 
of the full-sky Ce spectra. Secondly, it is clear from the simula- 
tions that the accuracy also depends on the cosmological models 
one wishes to investigate. If the parameter estimation procedure 
is restricted to a certain range of parameters (e.g. by ignoring 
cosmological models ruled out by previous experiments such as 
WMAP, 2dFGRS and/or SDSS), then the approximation may 
be used without too much loss of accuracy even for an experi- 
ment with a large sky-cut like QUaD. We are therefore confident 
that our approximation will be applicable to parameter estima- 
tion from forthcoming CMB temperature and polarization ex- 
periments with differing sky survey areas such as QUaD, BICEP 
and/or Planck. 

5 CONCLUSION 

Pseudo-Cf estimators and their use in reconstructing full-sky 
power spectra have become popular tools in recent years for 
measuring the temperature and polarization power spectra from 
mega-pixel observations of the CMB (Netterfield et al. 2002; 
Ruhl et al. 2002; Benoit et al. 2003; Hinshaw et al. 2003). In 
this paper, we have investigated in detail the reconstruction of 
the full set of CMB temperature and polarization full-sky power 



spectra from the pseudo-C, estimators. We have presented the 
calculation of the pseudo-C, method relating the observed to 
full-sky power and cross-power spectra for the complete set of 
CMB temperature and polarization spectra: Cj T , Cf E , C BB , 
Cf E , Cj B and Cf B . These calculations provide one with a 
mechanism to reconstruct all six CMB power spectra from ob- 
servations of the CMB over a finite region of sky. 

Having calculated the relevant coupling matrices required 
for the power spectra reconstructions, we have tested these re- 
constructions, which we have formulated in terms of bandpow- 
ers, on simulated CMB temperature and polarization data sets. 
We have conducted these tests on realistic simulations of two 
upcoming CMB experiments - (1) the Planck satellite where we 
have considered the reconstruction of the power spectra from 
two sky surveys over 14 months and (2) a combination of the 
ground-based QUaD polarization experiment with the temper- 
ature data from 4-year WMAP observations. In both cases, we 
have found the estimators to be unbiased over the entire range 
of multipoles investigated between £ = and £ — 2500. For 
these tests, we have included primordial gravitational waves and 
the effects of gravitational lensing in our input cosmological 
model in order to test the reconstruction of the B-mode polar- 
ization power spectrum. Using an input tensor-to-scalar ratio of 
T/S — 0.05, we recover a significant detection of the lensing- 
induced B-mode power from both the Planck and QUaD sim- 
ulations for multipoles, £ < 1000. The reconstructions of the 
non-cosmological cross-spectra, Cf B and Cf B are consistent 
with zero on all scales for both experimental scenarios. The re- 
construction of these two spectra can be useful as a test for sys- 
tematic effects in either the data itself or in the data analysis 
pipeline of CMB experiments. 

We have proposed a new method for fast and accurate cal- 
culation of the covariance matrix of the pseudo-C, estimators 
for any given cosmological model. Having presented the expres- 
sions for the full covariance of the pseudo-C, estimators, we 
have made a symmetrized approximation, similar to that used by 
Efstathiou (2004) for the temperature-only case. Our proposed 
method is based on fitting the mixing matrices of the approx- 
imated expressions from the covariances measured from sim- 
ulated data sets. We have tested this semi-analytic approxima- 
tion by comparing it with the covariance measured from simu- 
lations for a range of cosmological models drawn from within 
the WMAP 3(J parameter space region. These tests have been 
carried out for both a Planck-like experiment, with a galactic 
cut applied and for a QUaD-like experiment covering ~ 300 
deg 2 . We find that the semi-analytic approximation reproduces 
all the major features of the simulated covariances with varying 
degrees of accuracy which depend mainly on survey area and 
the range of cosmological parameters investigated. The effect of 
our approximation on the estimation of cosmological parameters 
will be investigated in a forthcoming paper. Comparing the er- 
rors measured from the simulations with those predicted by the 
"Knox" formulae, we find that the latter under-estimate the true 
errors by ~ 10 per cent for a near full-sky experiment and by 
~ 25 per cent for a small-scale experiment like QUaD which 
spans ~ 0.78 per cent of the sky. 

In the not-too-distant future, there is the prospect of high- 
quality full-sky mega-pixel maps of not only the CMB tempera- 
ture field, but also of the polarization field from satellite missions 
such as WMAP, Planck and later perhaps, NASA's proposed In- 
flation Probe experiment. More immediate prospects come from 
a number of novel ground-based polarization experiments such 
as QUaD and BICEP, both of which are due to begin observa- 
tions in 2005. Beyond this, the CLOVER experiment in planned 
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Figure 5. The coiTelation matrix of the pseudo-Cf measurements from the Planck-like simulations (lhs) and using the semi-analytic approximation 
(rhs). The input cosmological model is taken from the WMAP lcr parameter space region. For clarity, we plot the correlations for 1 ^ t ^ 100 only. 
Multipoles in the range, 1 ^ £ ^ 100 are the TT spectra, 101 < £ < 200 are EE, 201 < I ^ 300 are BB, 301 < £ < 400 are TE, 401 ^ £ < 500 
are TB and 501 ^ £ ^ 600 are the EB spectra. 
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Figure 6. The square root of the diagonal elements of the covariance matrices depicted in Fig. 5. In each plot, the light gray line is the error measured 
from the simulations, the red (dark gray) line is the semi-analytic approximation and the heavy dashed line is the error calculated using the "Knox" 
formulae. The errors plotted are (clockwise from top left) TT-TT, EE-EE, TE-TE, EB-EB, TB-TB and BB-BB. 
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Figure 7. The correlation matrix of the pseudo-CV measurements from the QUaD-like simulations (lhs) and using the semi-analytic approximation 
(rhs). The input cosmological model is taken from the WMAP la parameter space region. For clarity, we plot the correlations for 1 $C £ $C 100 only. 
Multipoles in the range, 1 < I < 100 are the TT spectra, 101 < £ ^ 200 are EE , 201 < £ ^ 300 are BB, 301 ^ £ ^ 400 are TE, 401 ^ £ 4, 500 
are TB and 501 ^ £ < 600 are the EB spectra. 
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Figure 8. The square root of the diagonal elements of the covariance matrices depicted in Fig. 7. In each plot, the light gray line is the error measured 
from the simulations, the red (dark gray) line is the semi-analytic approximation and the heavy dashed line is the error calculated using the "Knox" 
formulae. The errors plotted are (clockwise from top left) TT-TT, EE-EE, TE-TE, EB-EB, TB-TB and BB-BB. 
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Figure 9. The correlation matrix of the pseudo-C^ measurements from the QUaD-like simulations (lhs) and using the semi-analytic approximation 
(rhs). The input cosmological model is taken from the edge of the WMAP 2cr parameter space region. Multipoles in the range, 1 ^ I ^ 100 are the TT 
spectra, 101 < t ^ 200 are EE , 201 ^ £ ^ 300 are BB, 301 ^ I ^ 400 are TE, 401 ^ £ < 500 are TB and 501 ^ £ ^ 600 are the EB spectra. 
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Figure 10. The square root of the diagonal elements of the covariance matrices depicted in Fig. 9. In each plot, the light gray line is the error measured 
from the simulations, the red (dark gray) line is the semi-analytic approximation and the heavy dashed line is the error calculated using the "Knox" 
formulae. The errors plotted are (clockwise from top left) TT-TT, EE-EE, TE-TE, EB-EB, TB-TB and BB-BB. 
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Figure 11. The simulated (upper, left) and semi-analytic (upper, right) correlations for the WMAP lcr region Planck simulations shown in Figs. 5 and 
6. The residuals between the two, which have an rms value of ~ 0.03, are shown in the lower plot. Once again, for clarity, we plot the correlations for 
1 < I < 100 only. Multipoles in the range, 1 < i < 100 are the TT spectra, 101 ^ I s£ 200 are EE , 201 < I < 300 are BB, 301 < I < 400 are 
TE, 401 < I < 500 are TB and 501 < i < 600 are the EB spectra. 



for 2008. These present the astronomical community with the 
opportunity of continuing the rapid progress in cosmology al- 
ready achieved with the 2dFGRS, WMAP and SDSS programs. 
In addition to improving parameter constraints by up to a fac- 
tor of 6 (e.g. Bowden et al. 2004), such polarization information 
may also provide an opportunity to detect and perhaps measure 
the amount of primordial gravitational waves in the Universe 
through measurements of the B-mode polarization power spec- 
trum. Such a measurement would constitute an important addi- 
tion to our knowledge of the physics of early Universe, possibly 
providing a handle on the energy scale, and other parameters, 
of inflation (e.g. Kinney 1998). If the rapid progress in CMB 
observations continues as envisaged, then the methods and re- 
sults presented here can be used to quickly and accurately extract 
power spectra and cosmological parameters from CMB temper- 
ature and polarization data. 
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APPENDIX A: PSEUDO-C^ COUPLING MATRICES 



As discussed in Section 2.3, the pseudo-CV spectra measured from a finite region of sky are related to the underlying full-sky power 
spectra by the matrix relation, 
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for M,N = (0, +, — )• The W^,^ , matrices are defined as in equations (8), (11) and (12). In order to relate Ci to C e >, one needs to 
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numerically calculate expressions of the type 
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with sW^, m given by equation (8) or (12). These are cumbersome to evaluate due to the high number of summations over multipoles 
aggravated by the need to perform integrals over the whole sphere. Fortunately, the coupling matrix, M«/, can be simplified using 
the properties of the Wigner-3j symbols. One starts by expanding the window functions, Wt{Q) and Wp(£Y), in spin-0 spherical 
harmonics such that their auto- and cross-power spectra (which we have scaled here by 21 + 1) are given by 
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where the wf m and wf m coefficients are the spherical harmonic transform coefficients of the window functions: 

wj m = J daw T (n)Y; m , 

wL = J dQW P (Q)Y e * m . 
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In order to calculate terms like equation (A4), we apply the properties of integrals over solid angle of the spin-weighted spherical 
harmonics, 
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and finally use the orthogonality relation of the Wigner 3-j symbols (Varshalovich, Moskalev & Khersonskii 1998), 
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Then it is straightforward to notice that all the imaginary terms present in the coupling matrix, M«/, vanish because of the symmetries 
of the Wigner 3-j symbols, in particular 
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The imaginary terms are either associated with the expressions W° e , or W^,, (and their complex conjugates) in equation (A2). For the 
term W?,7 , we have a dependency of the type, 
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thus F(£, £',£") — 0. The same line of reasoning applies for the term W^T ■ Therefore, we obtain the simplified relation, 
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where S = T if s = and S = P if s = +2 or —2. The value of S' depends in the same way on the value of a'. Using the previous 
result and after some algebra, we obtain for the non-zero mode-mode coupling matrix M«/ components (see equation 19) the following 
set of expressions: 
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APPENDIX B: COVARIANCE MATRIX OF THE PSEUDO-ft ESTIMATORS 

In Section 3, we have described our procedure for calculating the semi-analytic approximation for the pseudo-CV covariance matrix. 
Here, we present the explicit expressions for the semi-analytic approximation. We write these in terms of the various mixing matrices, 
X^, cd , which we define by 
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where {a, b, c, d} can take any of the values {0, +, — } and the W matrices are the entries of the window matrix, equation (15), and are 
defined in equations (8), (11) and (12). In terms of these matrices then, the approximations to the various pseudo-Cf covariances are: 
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